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Abstract. We investigate f(R) theories of gravity within the Palatini approach and show how one can determine the expansion 
history, H(a), for an arbitrary choice of /(/?). As an example, we consider cosmological constraints on such theories arising 
from the supernova type la, large scale structure formation and cosmic microwave background observations. We find that best 
fit to the data is a non-null leading order correction to the Einstein gravity, but the current data exhibits no significant preference 
over the concordance ACDM model. Our results show that the often considered 1 /R models are not compatible with the data. 
The results demonstrate that the background expansion alone can act as a good discriminator between modified gravity models 
when multiple data sets are used. 



1. Introduction 

The combination of Einstein's General Relativity (GR) and 
ordinary matt er, as described by th e standard model of par- 
ticle physics ( lEidelman et alJl2004h . cannot explain the cur- 
rent cosmological observations. Key observations confronting 
the matter only universe are the luminosity-redshift rela- 
tionship from observations of supernovae of type la (SNIa) 
(j^iess et al. 2004), the matter power spectrum of large scale 
structure as inferred from gala xy redshift surveys like the 
Sloan Digital Sky Survey (SDSS) ('Xegmark et al."2004') and the 
2dF Galaxy Redshift Survey (2dFGRS) (Colless et al. 2001), 
and the anisotropie s in the Cosmic Mi crowave Background 
Radiation (CMBR) fe pergel et al^ l2003^. In order to account 
for the results from all of these cosmological probes within 
GR, two exotic components are required in the matter-energy 
budget of the Universe. These two components are respectively 
dark matter, a collisionless and pressureless fluid, which con- 
tributes about 25% of the universe energy budget, and a neg- 
ative pressure fluid called dark energy. Currently the dark en- 
ergy component dominates the energy density of the universe, 
causing accelerating expansion. 

Despite the very good agreement between the so-called 



one can ask whether the necessity of including dark matter and 
dark energy in the energy budget is a sign of our lack of un- 
dersta nding of gravitational physics (Lue et al. 2004a- b.: iLud 
l2003l) . A natural alternative to adding new exotic fluids is to 
modify gravitational physics. Several wa ys of mo difying g rav- 
ity in order to dispense with dark matter jSanders & McGaughl 
2003 iMilgromI fT994l iMoffatI l2004 iBekensteirJ 1200.4 
Skordi s et al.' "2006; 'Sellwood & Kosowsky' 7001) or dark 

urneJ 



enerj 
5003 



rv CDvali et aL .2000.: Deffavet 2001: Dvali & 
lArkani-Hamed et alJ l2002t iNoiiri & Odintso' 



'2004aV, '2005"; 'Abdafla et al.' '2005'; "Freese & Lewis' "2002; 
Meng & Wang 2003.: Dolgov & Kawasaki 2003: Shao et alj 




can 
vations 



&J 12004^ lAh med et alJ l2004 

have been proposed, and they 
with some de gree of success account for the obser- 

, JT.ue et alJ ^04^ ^ Eooj iBento et alJ Eooj 

"2003; Deff^ayet et al 2002; Carroll et al. 2005; Multamaki et al.l 
. 2003. 2004: Koivisto et al, 2005: Elgarov & Multamaki 2005: 
lAmarzguioui et al j2005U]VIultamaki & Elgaroyl2004 . 

In this article we investigate a family of alterna- 
tive models to the dark energy paradigm based on 
a generalization of the Einstein-Hilbert Lagrangian. 
These models are call ed Nonli near Theories of Gravity 



dark matter and dark energy is one of the greatest myster- 
ies of modern cosmology. In fact, none of the dark mat- 


theories (Bronnikov & Chernakova 2005: Cognolaet alJl2005l 


Nunez & Solganik 


2004 




Ezawa et al. 


2003 




Barraco et al.' 



ter candidates from high energy physics beyond the standard 


model (BertoneetalJl2005: Ellis 2000; 


Brookfield et al. 


2005; 



lAmendo la 2000) have ever been observed. One should bear in 
mind that the existence of the dark sector is only inferred from 
the motion of ordinary matter in a gravitational field. Hence 



2000; Schmidt .1998: JU pplet aL. .1996: .BaiTow & Ottewilli 
1983), since the scalar curvature R in the Einstein-Hilbert 
Lagrangian is replaced by a general function f(R). The main 
motivation for this generalization is the fact that higher order 
terms in curvature invariants (such as R^, R^^R^y, R''^"^R^vais, 
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etc) have to be added to the effective Lagrangian of the 
gravitational field when quantum corrections are considered 
(iBuchbinder et al.' '1992'; 'Birrell & Davies| Il982t IVilkoviskvl 
[l992: Gasp erini & Veneziano 1992). Furthermore, there is no 
a priori reason to restrict oneself to the simple Einstein-Hilbert 
action when a more general formulation is allowed. Higher 
order terms in the gravitational action also have interesting 
consequences i n cosmology, hke riatural inflationarv behavior 
at early times ('Starobinsky"l980'; ^ Barrow & Cotsa ki'i '1988', 
Il99lt ILa & Steinhardt 1989), a nd late time acceleration 
of th e Universe (Meng & Wan^ l2003l l2004t ICarroU et alJ 
1200-4 ICanozziello et aL 120051) . Consequently, several au- 
thors have investigated whether such theories are indeed 
compatible with current cosmological observations, big 
bang nucleosynthesis, and solar system constraints in its 
weak field hmit (Clifton & Barrow 2005; Barrow & Clifton 
2006| fOuandt & Schmidt 1991; Dominauez & Barraco 
2004^ "Sotiriou 2005a b; Olmo 2005a.b: Canozziello 2002; 
parloni et al. 2005; Barraco et alJl2002t iHw ang & Noh 2001; 
Mena et alJ 120051: lAllemandi et alJ l2004albt l2005l) . Most of 
these investigations have been model dependent and the 
conclusions are somewhat contradictory (Flanagan 2004; 
IVollickl2004h. 

When dealing with f{R) theories of gravity, the choice of 
the independent fields to vary in the action is a fundamental 
issue. In the so-called Palatini approach one considers the met- 
ric and the connection to be independent of each other, and 
the resulting field equations are in general different from those 
one gets from varying only the metric, the so-called metric ap- 
proach. The two approaches lead to the same equations only if 
f(R) is linear in R. The correct choice of approach to derive the 
field equations is still a hot topic of research. Initially, f{R) the- 
ories were investigated in the metric approach. However, since 
this method leads to fourth-order equations and the Palatini ap- 
proach leads to second-order equations, the latter is appeal- 
ing because of its simplicity. Moreover, the equations result- 
ing from the metric approach see m to have instability prob- 
lems in ma ny interesting cases dPolgov & Kawasaki! l2003l 
IChiball200/!h from which the Palatini approach does not suf- 
fer. However, recent work has cast do ubts over these instabili- 
ties JCembran os 2005 : Sotiriou 2005a!). In the present paper we 
will concentrate on the Palatini approach. 

The aim of this article is to use current cosmological 
data to consider possible deviations from GR by combining a 
number of different cosmological observations. The data used 
are the latest Supernov ae la gold s et (Riessetal. 2004), the 
CMBR shift parameter (Bond et al. 1997), the baryon oscilla- 
tion length scale cEisenstein et aL.2005.) and the linear growth 
factor at the 2dFGRS e ffective redshift {Hawkins et al. 20031 
IWang & Teg mark 2004). As far as we are aware, this is the first 
time one uses all of the main cosmological data sets in order to 
constrain these models. 

The structure of the paper is as follows: In section |2] we 
deduce and summarize the basic equations and properties of 
general f{R) gravities in the Palatini approach. In section|3l we 
investigate observational constraints based on the evolution of 
the background of the Universe. In particular, we consider fits 
to the SNIa and CMB shift parameter from the WMAP data. In 



section 13 we analyze the evolution of linear perturbations and 
probe large scale structure formation in these models using the 
linear growth factor derived from the 2dFGRS data. Finally, 
section|5lcontains a summary of our work and our conclusions. 

2. General f{R) Gravity theories 

The action that defines f{R) gravity theories in the Palatini for- 
malism is the following: 



(1) 



where k - SttG, S mlgfiv^^f'm] is the matter action which de- 
pends only on the metric gfjy and on the matter fields t//,„, 
R = R{g, F) - g^^Rfjyit) is the generalized Ricci scalar and 
Rftv(r) is the Ricci tensor of the affine connection F, which in 
the Palatini approach is independent of the metric. The gener- 
ahzed Riemann tensor is given by ( Vollick 2003): 



pA pa 



(2) 



We define the Ricci tensor by contracting the first and the third 
indices of the Riemann tensor. The field equations will be ob- 
tained using the Palatini formalism i.e. we vary both with re- 
spect to the metric and to the connection. 

Varying the above action with respect to the metric, we ob- 
tain the generalized Einstein equations 



(3) 



where f'(R) = df/dR and is the energy momentum tensor 

2 6S 



T = - 



(4) 



Varying with respect to the connection f and contracting gives 
us the equati on that determines the generalized connection 
llVonickl200i : 



^a^f'(R)^F88n-0, 



(5) 



where V is the covariant derivative with respect to the affine 
connection F. This equation implies that we can write the affine 
connection as the Levi-Civita connection of a new metric /i^v = 
f'(R)gf,v The Levi-Civita connections of the metrics g^jy and 
hjjy are then related by a conformal transformation. This allows 
us to write the affine connection as 



(6) 



where F^^ is the Levi-Civita connection of the metric g/jy. 
The generalized Ricci tensor can now be written as 

„ 3 V,/-Vv/' , V,Vv/' , 1 WV.f 
V = l^f^vig) - 2 — ~2 + f, + 2^^" // ■ 

Here Rf,y(g) is the Ricci tensor associated with gfjy, and the 
covariant derivative associated with the Levi-Civita connection 
of the metric. 
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Since we are interested in cosmological solutions we con- 
sider the spatially flat FRW metric, 

ds^ = -dt^ + a{tf-5ijdx'dxK (8) 

and the perfect fluid energy-momentum tensor — 
diag(-p, p,p,p). 

Taking the trace of Eq. Q, gives 

Rf'(R)-2f(R)^-KT, (9) 

where T = g'^^ T^v - -(p - 3p). Eq. ^ can in certain special 
cases be solved explicitly for R - R(T), but this is not possible 
in general. 

The generalized Friedmann equation can be derived 
straightforwardly using the generalized Ricci tensor, Q: 



H + 



If'Y lK{p + 3p) If 



2f'j 6 /' 6f 



(10) 



which a grees with the result found in lAflemandi et alJ 
b) (note the sign difl'erence in the definition of the 
Einstein equation)'. 

If the equation of state of the fluid, p = pip), is known, 
one can use the continuity equation of the fluid together with 
Eq. to express R as 



R = -3H 



{l-3p'(p))(p + p(p)) 
Rf"(R)-f'(R) 



(11) 



Using the Friedmann equation it is easy to see that now H = 
H(p, R), which together with Eq. forms an algebraic set of 
equations from which one can in principle determine H - H{R) 
for any fluid. 

In particular, in the case of a constant equation of state, p = 
wp, we get Kp = (R f - 2/)/(l - 3w) and hence the generalized 
Friedmann equation can be written as 

^^2^ 1 (I + 3yv)Rf'- 3(1 +w)f 
-6(l-3w)/'(^_3^^,^^/W_^)- 

Similarly, from Eq. (|9j, we can write a - a{R): 



1 \ iU^) 

iRf'-2f)] 



Kpoil - 3w) 



(13) 



where po = p(ao) and we have chosen aq = 1 ■ Using these two 
equations, J12l i and il3i . it is now straightforward to determine 
the expansion history, H{a), for any f{R). 

In the special case w - 0, applicable when the universe is 
matter dominated, these equations reduce to 



1 Rf - 3/ 

6/' /, _ 3 /"(«/' -2/) f 

V 2 /'(«/"-/'>; 



a ^ { — (Rf'-lf)] . 



(14) 



(15) 



This differs from the results of lWang & Mend J2004 bv 3Hf'/f'. 



2.1. The leading correction to Einstein gravity 

In order to investigate to what extent observations allow devia- 
tions from general relativity (where f(R) = R), we consider the 
following gravity Lagrangian: 



f(R) = R 



1 + a 



R 



Hi) 



(16) 



where a, j3 are dimensionless parameters (note that in out nota- 
tion R is negative). Specializing to a matter dominated universe 
(w - 0), we get from Eq. (|9jl the relation between the curvature 
scalar and matter density: 



Kpm = -R 



l+ai2-/3) 



R 



Hi/ 



(17) 



We wish to recover standard behaviour at early times. This im- 
plies that the correction term in the Lagrangian must vanish for 
large \R\, and hence we must demand that ^0 < 1. Furthermore, 
we must demand that the right hand sides of equations il4\ 
and in\ always remain positive, which restricts the parameter 
space further. 

Defining Q.,„ = /fpJJ,/(3//,)), and choosing units so that Hq - 
1 , we can solve for Rq from Eq. ( I17> . Consistency then requires 
that substituting the obtained value of Rq into Eq. il2\ must 
give //() = 1 . Hence, given a and /3, Q.„, is fixed. 

As an example, consider the case y6 = 0, which corresponds 
to flie ACDM model. From Eq. O, we have Hi = (3aHl - 
2Ro)/6 and from Eq. J17l i. 30,„//y - 2aHl - Rq, and hence 

- I + a/6 is fixed. 

3. Observational constraints from the Background 
Evolution 

Armed with the modified Friedmann equation, we can now 
consider the constraints arising from cosmological observa- 
tions. In this section we will consider quantities related to the 
background expansion of the Universe: the SNIa luminosity 
distance-redshift relationship, the CMBR shift parameter and 
the baryon oscillation length scale. 

3.1. CMBR Shift Parameter 



The CMBR shift parameter ( Bond et al. ll997':'M elchiorri et alJ 
■Odman et al..2003) in a spatially flat universe is given by 



^Q.,nH^ 



dz 



(18) 



where Zder is the re dshift at decoupling. The WMAP team 
dSpergel et alJ l2Q03) quotes Zdec = 1088!' and K = 1.716 ± 
0.062. In writing the shift parameter in this form we have im- 
plicitly assumed that photons follow geodes ies determined by 
the Levi-Civita connection. In Koivisto (2005) this is shown to 
be the case if there is no torsion present. Furthermore, in order 
to use the shift parameter, the evolution of the universe needs 
to be standard up to very late times so that at decoupling we re- 
cover standard matter dominated behaviour Hence, we restrict 
our analysis to/3 < I. 
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Fig. 1. The 68, 95 and 99% confidence contours arising from 
fitting the CMBR shift parameter. The parameter values corre- 
sponding to the concordance ACDM model (Q,„ = 0.27, Qa = 
0.73 or Of = -4.38, y6 = 0) are marked with a cross. The gray 
area represents a section of the parameter space that is not al- 
lowed. 




-10 -9 -a -7 -6 -S -4 -3 -2 -1 
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Fig. 2. The 68, 95 and 99% confidence contours arising from 
fitting the SN la data. The parameter values corresponding to 
the concordance ACDM model (Q,„ = 0.27, Qa = 0.73 or 
a = -4.38, y6 = 0) are marked with a cross. The gray area 
represents a section of the parameter space that is not allowed. 



Since we do not have an explicit expression for the Hubble 
parameter in terms of the redshift, it is useful to rewrite the shift 
parameter in terms of the curvature scalar, R: 

dz 



I — 2 r*" 



Hiz) 
(R) dR 



(R)2 H{R) 

**' Rf" - f 



dR 



{Rf'-2ffl^ HiR) 



(19) 



The constraints arising from the CMBR shift parameter can 
be seen in figure [2 




3.2. SNIa constraints 

In order to incorporate measurements from SNIa, it is useful to 
rewrite the expression for the luminosity distance as 

dz 



Ddz) ^ (l+z) 



H(z) 



(20) 



- ^-^lnJiiiRf'-2fy''x 
Rf" - f dR 

J«„ {Rf'-2ff/' HiRY 

For the supern ova data, we use the "Gold data set" from 
iRiess et a T (2004). The contour plots showing the constraints 
on the parameters from the supernovae can be seen in figure|3 
To get these plots we marginalized over the Hubble parameter 
h. 

With the added information from the CMBR in the form of 
the shift parameter, the situation improves as one can see from 



Fig. 3. The combined 68, 95 and 99% confidence contours 
arising from fitting the SN la and the CMBR shift parameter 
data. The parameter values corresponding to the concordance 
ACDM model (Q,„ = 0.27, Qa = 0.73 or or = -4.38, yS = 0) 
are marked with a cross. The gray area represents a section of 
the parameter space that is not allowed. 

from figure|3 Still quite a large degeneracy persists on the 99% 
level, but on the 68% level, the model is quite well constrained 
and centered around the concordance ACDM model. 



3.3. Baryon Oscillations 

The baryon oscillations in the galaxy power spectrum are im- 
prints from acoustic oscillations prior to recombination, which 
are also responsible for the the acoustic peaks seen in the 
CMBR temperature power spectrum. The physical length scale 
associated with the oscillations is set by the sound horizon at 
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recombination, whic h can be estimated from the CMBR data 
dSpergel et alJ2003h . Measuring the apparent size of the oscil- 
lations in a galaxy survey allows one to measure the angular 
diameter distance at the survey redshift. Together with the an- 
gular size of the CMB sound horizon, the baryon oscillation 
size is a powerful probe of the properties and evolution of the 
universe. 

The imprint of the primordial baryon-photon acoustic os- 
cillations in the matter power spectrum provides us there- 
fore wit h a 'standard ruler' via the dimensi onless quantity A 
llLindd 120031 boosi IH u & Sugi vamalll994 lEisenstein & Hul 
ll998HEisenstein & WhitQ.20041) : 



jQ.mE(zi) 



-1/3 



1 r Ji- 

z\ Jo E{z) 



(21) 



where E(z) = i/(z)/i/o. 

Recently the acoustic signature associated with the bary- 
onic oscillations has been identified at low redshifts in the dis- 
tribution s of Luminous Red G alaxies in the Sloan Digital Sky 
Survey (lEisenstein et alJ2005l) . with a value of 



where 



A = D,(z = 0.35) 



0.35c 



0.469 ±0.017, 



cz 



Hiz) 



1/3 



(22) 



(23) 



and Dm(z) is the comoving angular diameter distance. For in- 
stance, in the case of ACDM with Q,„ = 0.3, Qa = 0.7 and 
h = 0.7 we have D,(z = 0.35) = 1334Mpc. 

Using the observed baryon oscillation length scale, one can 
hence constrain the cosmological model. The confidence con- 
tours for the modified gravity model we are concerned with in 
this paper, are shown in figure^] 

Once again there is a large degeneracy in the a vs /3 plane, 
similar to the case when we fit SNIa (figure |2j and the CMBR 
shift parameter (figure [Q. Such degeneracies are strongly re- 
stricted, however, when one combines all of the data sets in 
one single plot (see figure|5}. 

Using the baryon oscillations length scales and combining 
it with the SNIa and the CMBR shift parameter data hence im- 
poses strong constraints on possible deviations from Einstein's 
General Relativity. This also demonstrates how combining the 
current data, one can efficiently study and constrain cosmolog- 
ical models from the background expansion only. 

The best fit model to the three data sets is a = -3.6 and 
y6 = 0.09. This is sUghdy different from the ACDM model, but 
the ACDM model is well within the 68% confidence contour. 



4. Large Scale Structure: Formalism and 
Constraints 

So far we have only considered observables related to the back- 
ground evolution. In order to gain further information, it is use- 
ful to go beyond these 'zeroth order' tests and consider per- 
turbations within f{R) models. Knowledge on the evolution of 




Fig. 4. The 68, 95 and 99% confidence contours arising from 
fitting the length scale associated with the baryon oscilla- 
tions. The parameter values corresponding to the concordance 
ACDM model (Q,„ = 0.27, Qa = 0.73 or o- = -4.38, yS = 0) 
are marked with a cross. The gray area represents a section of 
the parameter space that is not allowed. 




Fig. 5. The 68, 95 and 99% confidence contours arising from 
fitting the combined data from the SNIa, CMBR shift pa- 
rameter and the length scale associated to the baryon oscilla- 
tions. The parameter values corresponding to the concordance 
ACDM model (Q„, = 0.27, Q.^ = 0.73 or a = -4.38, yS = 0) 
are marked with a cross. The gray area represents a section of 
the parameter space that is not allowed. 



perturbations allows one to confront models with large scale 
structure observations from galaxy surveys. 

Several authors have investigated cosmological pertur- 
bations in gen eralized gravity theories in the metric ap- 
proach ( iHwangt .1991 a.b>) . However within the Palatini for- 
malism, only very recently the first steps were made 
( iKoivisto & Kurki-Suoni o 2005). Here we wi ll follow a spher- 
ical collapse formalism (Lue et al.. .2004at iMultamaki et alJ 
.2003.) where by requiring that a general gravity theory respects 
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the Jebsen-Birkhoff theorem^ one derives the modified grav- 
itational force law necessary to describe the evolution of the 
density perturbations at astrophysical scales. 

4.1. Spherical collapse in f(R) theories 

Although the spherical collapse model has been set and used 
for a long time (Peebles 1993; Padmanabhan 1993) and in 
many d i fferent contexts dLue e t al. 2004a; Mota & B arrow! 
l2004albt iMota & van de Bruckl l2004t iClifton et al...200 5'). it 
ha s not been applied previously to Nonlinear Gravity Theories. 
In lLue et alJ (l2004ah . one starts by assuming a generalization 
of the Jebsen-Birkhoff theorem: for any test particle outside 
a spherically symmetric matter source, the metric observed 
by that test particle is equivalent to that of a point source of 
the same mass located at the center of the sphere. With this 
one assumption, one can deduce the Schwarzschild-like met- 
ric of the new hypothetical gravity theory. Armed with the 
Schwarzschild-like metric one can then investigate the evolu- 
tion of spherical matter overdensities and compare it to the 
latest large scale structure data. The idea is the following: 
Consider a uniform sphere of dust. Imagine that the evolu- 
tion inside the sphere is exactly cosmological, while outside 
the sphere is empty space, whose metric (given the Jebsen- 
Birkhoff theorem) is Schwarzschild-like (as defined by the met- 
ric equation (2.2) in Lue et al. (2004a)). The mass of the matter 
source (as determined by the form of the metric at short dis- 
tances) is unchanged throughout its time-evolution. The surface 
of the spherical mass therefore charts out the metric through all 
of space as the sphere expands with time, as long as we de- 
mand that the cosmological metric just inside the surface of 
the sphere smoothly matches the Schwarzschild solution just 
outside. In order to see how the metric depends on the mass 
of the central source, we just take a sphere of dust of a differ- 
ent initial size, and watch its surface chart out a new metric. 
The procedure for determin ing the metric from the cosmologi- 
cal evolution is described in iLue etalJ ( E004al) and we refer the 
reader to this article for further information and details. 

An open question is the validity of the Jebsen-Birkhoff the- 
orem in f(R) theories. Although one can not explicitly show 
that the only possible solution to the field equations, when 
a spherically symmetric ansatz is inserted into them, is the 
Schwarzschild metric, there is strong evidence that this is in- 
deed the case. In fact that was show n for the case of f{R) = 
R + i f-, where R^ = R'^ a'^'R" Pfiv JRamaswamv & Yasskirj 
11979 *). And it was generalized later for any type of invariant 
of the form R^, even in the case of a non-null torsion (Neville 
The authors also claimed that similar results would most 
probably be valid even for the case of higher order curvature 
invariants, such as R^, R*, etc. However, there is no math- 
ematical proof of this as of yet, even though there are sev- 
eral studies and proofs for other complex cases such as multi- 
dimensional gravities, Einstein- Yang-Mills systems, and con- 



^ What is commonly known as BirkhofF's theorem was first probed 
by the Norwegian J. T. Jebsen in 1921, two years prior to BirkhofF's 
work. See_physics/0508 1 63 /or more historical details. 



formally transformed metrics ( iBrodbeck & StraumannI 1 199*31 
iBronnikov & Melnikovl 19951) . 



4.1 .1 . Growth of perturbations 

We want to follow the evolution of a lump of matter throughout 
the history of our Universe. We start by considering a Universe 
with a Schwarzschild-like metric. 



ds^ = goo(f, r)dT^ - grrit, r)dr^ - r^dQ. , 



(24) 



where the metric componen ts are uniquely determined from a 
given a(t) jLue et alJ2004ah : 



goo = £'[1 - i-fa^] 
a 



(25) 



where E - sjgoogrr is a constant. 

In order to investigate the evolution of density perturbations 
in thi s scenario we use the spherical collapse model ( Peebles, 
Il993h . Consider a top-hat overdensity 6(t) of a spherical distri- 
bution of dust with mass M and radius r defined by 



1+6 



M 



(26) 



where p(f) is the background matter density. Using the geodesic 
equation as expressed by differentiating Eq. i25i with respect 
to f, we get 



Id, r, 



(27) 



where g{x) is defined via a generalized Friedman equation 

^(^) = |J, (28) 

and where we have expressed p in terms of a dimensionless 
quantity x defined as 



X = 



(29) 



Using the constraint from the Jebsen-Birkhoff theorem, one 
can calculate the evolution of an overdensity by following the 
geodesic of a spherical mass, without the need to consider what 
is happening outside the spherical mass itself. This is only pos- 
sible if spherically symmetric configurations respect the metric 
equations i25\ (Lue et al. 2004a). Differentiating equation i26\ 
twice with respect to time, and using equation i21\ . one obtains 
a new equation for (5(f) 

.41., 
6 + 2H6--- 



3(1 +6)Hl 



31+6 

^x(l+6)g'(x(l+6))-g(x(l+6)) 



-3(1 + 6)Hl 



-xg'(x) - g(x) 



(30) 
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Fig. 6. The 68, 95 and 99% confidence contours arising from 
fitting the Hnear growth factor using the spherical collapse 
approach. The parameter values corresponding to the con- 
cordance ACDM model (Q.,„ = 0.27, = 0.73 or a = 
-4.38, y6 = 0) are marked with a cross. The gray area repre- 
sents a section of the parameter space that is not allowed. 



At linear order in perturbation theory, equation i30\ gives: 

6 + 2H6 = 4-nGp6[g'(x) + 3xg"(x)] . (31) 

In our case both H and p are expressed as functions of the scalar 
curvature R and not the time explicitly. We therefore need to 
rewrite (13 1> in terms of derivatives of R. The equation to be 
solved will then be of the form 



d^6 d6 
— ^ +AiR) — 
dR^ dR 



B(R)6 . 



(32) 



where A(R) and B(R) are two rather unattractive functions of 
the scalar curvature. Using Eq. i3H . one can solve for the evo- 
lution of the linear growth factor and compare to observations. 

4.1 .2. Constraints 

The large scale structure information we choose to use here is 
the linear growth rate F(zmf) = 0.51 ± 0. 1 1 measured by the 
2dFG RS iverde et alJ [2002: Kn op et alJ l2003: Hawkins etaD 
where F = d\nD/dlna. We compare the theoretical 
value we get for the linear growth rate for our model with 
the value measured by the 2dFGRS at its effective redshift, 
Z2dF = 0.15. The constraints arising from the linear growth 
rate are plotted in figure |6l Combining this with all the other 
constraints leads to the confidence contours shown in Fig.0 

5. Conclusions 

We have investigated observational constraints on f{R) theo- 
ries within the Palatini formalism. In order to relate these the- 
ories to observations, we have shown how one can determine 
the expansion history for a given /(R). In particular, in a mat- 
ter dominated universe, determining H{a) is straightforward as 
expressed by Eqs il4\ and il5\ . 



Fig. 7. The 68, 95 and 98% confidence contours arising from 
fitting to the Combined SNIa, shift parameter, baryon oscil- 
lations and the linear growth data sets. The parameter val- 
ues corresponding to the concordance ACDM model (Q^ - 
0.27, Qa = 0.73 or Q- = -4.38, /? = 0) are marked with a cross. 
The gray area represents a section of the parameter space that 
is not allowed. 



In particular we have investigated the possible form of 
the leading correction to standard GR, parameterizing the 
Lagrangian for gravity as - R- a^-RjH^^, and have used 
a combination of data sets to determine the allowed ranges of 
a and y6. This is by no means an exhaustive study. Other inter- 
esting forms of f{R) that are definitely worth studying include 
f{R) ^ \n (R) and perhaps especially f{R) ^ R-ci/R + czR^ 
(ISotirioiJ2 005b )■ but our main purpose here has been to set up 
the formalism and demonstrate the effectiveness of combining 
the current data sets. 

Using a combination of data sets that probe the background 
evolution, we have found that the cuiTent data efficiently con- 
strains the allowed parameter space of the leading correction 
to GR in the Palatini approach. The best fit models to the in- 
dividual data sets are (a,/3) - (-10.0,-0.51) for the SNIa, 
(a,/3) = (-8.4,-0.27) for the CMBR shift parameter and 
(a,/3) - (-1.1,0.57) for the baryon oscillations. The best fit 
to the combination of these data sets is {a,/3) - (-3.6,0.09), 
but the ACDM concordance model is well within the Icr con- 
tour Note, however, that the commonly considered 1 /R model 
is strongly disfavoured by the data. 

In order to bring in additional information from the cur- 
rent galaxy surveys, we have also considered the growth of 
structures in these models of modified gravity. By assuming 
that the new gravitational physics obeys a limited version of 
the Jebsen-BirkhofF theorem, we can describe the evolution of 
overdensities in f(R) gravity theories at sub-horizon scales. We 
find the best fit model to the linear growth factor alone to be 
ia,/3) - (-4.25,0.05) but the allowed parameter range is de- 
generate and does not improve constraints derived from the 
background evolution. In order to fully utilize the information 
available from the galaxy survey in form of the large scale mat- 
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ter power spectrum, a more detailed analysis is need ed along 
the lines presented in jKoivisto & Kurki-Suoniol2005l) . 

In summary, modified gravities provide us with an interest- 
ing alternative to the cosmological concordance model with a 
dominant dark energy component. Modern cosmological data 
can efficiently constrain such models. These data indicate that 
currently there is no compelUng evidence for non-standard 
gravity. 
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